o 

(N 



< 
oo 



o 
in 

o 



c3 



Extended Object Tracking with 
Random Hypersurface Models 

Marcus Baum and Uwe D. Hanebeck 

Intelligent Sensor-Actuator-Systems Laboratory (ISAS) 

Institute for Anthropomatics 

Karlsruhe Institute of Technology (KIT), Germany 

marcus . baum@kit . edu , uwe . hanebeck@ieee . org 



Abstract 



The Random Hypersurface Model (RHM) is introduced that allows for estimating a shape approx- 
imation of an extended object in addition to its kinematic state. An RHM represents the spatial 



extent by means of randomly scaled versions of the shape boundary. In doing so, the shape param- 
eters and the measurements are related via a measurement equation that serves as the basis for a 
Gaussian state estimator. Specific estimators are derived for elliptic and star-convex shapes. 



^ 1. Introduction 

In target tracking applications [lj where the resolution of the sensor device is higher than 
the spatial extent of a target object, the usual point object assumption is not justified as several 
,— i- different points, i.e., measurement sources, on the target object may be resolved during a single 

scan (see Fig. IT]). The resolved measurement sources typically vary from scan to scan and their 
locations depend on the shape of the object but also on further properties such as the surface or 
the target-to-sensor geometry. In this article, the basic idea is to approximate an extended object 
with a geometric shape such as an ellipse [21 El H] as depicted in Fig. [I] The tracking problem 
then consists of estimating the shape parameters in addition to the kinematic parameters. The 
locations of the measurement sources are not explicitly estimated. Reasonably, the shape of the 
target should be described as detailed as possible. However, when the measurement noise is rather 
high and only a few measurements are available, it may only be possible to infer a coarse shape 
approximation such as a circle. 

The unknown locations of the measurement sources are usually modeled with a probability 
distribution whose mass is concentrated on the extended object (a so-called spatial distribution) 
[21 Ej - In general, no closed- form solutions for the likelihood function resulting from a spatial 
distribution model exist, so that Monte Carlo methods are frequently used for approximating the 
Bayesian filter solution [21 El El El El 121 EH] ■ In case an elliptic extent is represented with a covariance 
matrix [3l llll[12lll3|, ll4 t ll5|. closed-form expressions can be derived with the help of random matrix 
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X = Measurement 



Figure 1: Shape approximation of an extended object with an ellipse. 

theory. Spatial distributions have been embedded into the Probability Hypothesis Density (PHD) 
filter for tracking multiple extended objects [161 El OS IIS1 EQ]. In [211 E2] it is proposed to drop 
all statistical assumptions on the measurement sources, which leads to combined set-theoretic and 
stochastic estimator. 

1.1. Contributions 

The main contribution is a novel systematic approach for modeling the unknown location of 
a measurement source on a spatially extended object called Random Hypersurface Model (RHM). 
The basic idea is to assume the measurement source to lie on a scaled version of the shape boundary, 
where the scaling factor is modeled as a random variable. In this manner, it is possible to form a 
suitable measurement equation that serves as the basis for constructing a Gaussian state estimator. 
In order to illustrate the novel approach, specific RHMs and corresponding Gaussian estimators 
are developed for ellipses and free-form star-convex shapes. To the best of our knowledge, this is 
the first extended object tracking method for explicitly estimating a free-form star-convex shape 
approximation. Actually, with this method, it is possible to track a target whose shape is a priori 
unknown and estimated from scratch over time. 

Remark 1 . This article is based on [231 [2H [251 IH ISH [271 I2S1 [291 ETi] ■ 

1.2. Overview 

The remainder of this article is structured as follows: In the following section, the general 
probabilistic framework for extended object tracking is presented. The new target extent model 
called Random Hypersurface Model (RHM) is subsequently introduced in Section pi Based on 
these models, a formal Bayes filter for extended object tracking is described in Section [4j Then, 
particular implementations of RHMs for ellipses (Section [5]) and star-convex shapes (Section [6]) are 
developed. Both shape representations are evaluated by means of typical extended object tracking 
scenarios in Section [3 This article is concluded in Section [8) 

2. Modeling Extended Targets 

The state vector of the extended object at discrete time k is represented with a random vectoiQ 
■HLk = [ m ft ' G*ifc) T >pT] that consists of the target location m k , a shape parameter vector p , and 
an optional vector x^, for the kinematics, e.g., velocity. The shape parameter p specifies a two- 
dimensional set that is denoted with S(p ) C K 2 . For example, a circular shape can be specified by 



1 In this article vectors are underlined, e.g., a; is a vector, and random variables are written in bold, e.g., x denotes 
a random vector. 
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Figure 2: Measurement model for a single measurement. 

"k, i.e., p" = r^ and the corresponding shape is S(p k ) = {z \ z G R 2 and \\z_W2 < r fc}n 
Note that we focus on two-dimensional shapes. Nevertheless the concepts are also applicable to 
higher dimensional shapes in the same manner, see for example |31j . 

2.1. Measurement Model 

At each time step k, a set of nt two-dimensional point measurements {y / }"= 1 of the extended 
object is available. We assume the measurements to be mutually independent for given target state; 
hence, a measurement model for a single measurement is sufficient. 

Extent Model. For a given state x k , the target extent model specifies a single measurement source 
z.k,i ^ rn k + ^(Pfc) on ^ ne extended object (see Fig. pi PI In the following Section pi we introduce 
the so-called Random Hypersurface Model for the target extent. 

Sensor Model. For a given measurement source z kl , the sensor model specifies the measurement 
y . We focus on Cartesian point measurements corrupted with additive Gaussian noise according 
to' 

y kl = z.k,i + Vk,i » (!) 

where the noise term v_ k t is zero- mean white Gaussian noise with covariance matrix T, k , . 

2.2. Dynamic Model 

In contrast to a point target, the temporal evolution of both the shape and kinematic parameters 
has to be modeled for an extended object. In this article, we focus on linear motion models according 
to 

x.k+1 = A kX k + w k , (2) 

where A^ is the system matrix and w k is white Gaussian system noise. 

3. Random Hypersurface Models 

In this section, a new target extent model called Random Hypersurface Model (RHM) is intro- 
duced, which describes the location of a single measurement source in a star-convex regionj^J 

Definition 1 (Star-Convexity). A set S C IR 2 is star-convex with respect to the origin iff [0, 0] £ 
S and each line segment from [0, 0] to any point in S is fully contained in S. 



2 The ci in p" indicates that it parameterizes a circle. 
3 Set operations "+" and "•" are meant element-wise. 



4 The concept of an RHM is much more general than presented here. The restriction to (two-dimensional) star- 
convex shapes is done for simplifying the following discussions. 
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Figure 3: Random Hypersurface Model (RHM) for an ellipse. 

3.1. Motivation: Implicit Measurement Equation 

The objective is to form an equation that relates the measurement source with the shape 
parameters. When the measurement source z k t lies on the boundary of the object, this is easy as 
the boundary is a closed curve (in general a hypersurface) that can be described by an implicit 
equation in the form g(z k i,p,) = 0. However, the measurement sources may also lie in the interior 
of boundary. In order to cover the interior, the basic idea is to scale the object boundary as 
described in Fig. [3| As the corresponding scaling factor s kj i & [0, 1] for a measurement source z kl 
is unknown, we model it as a random variable and treat it as an additional noise term. By this 
means, we obtain the implicit relation 

9(z k ,i,P k ,s k ,i) = (3) 

which forms, together with (fil), an implicit measurement equation (note that g(z k i,p ,1) = 
g(z k i,p )). In this sense, modeling a two-dimensional region is reduced to modeling a curve 
by means of the scaling factor. 

3.2. Definition 

According to the above motivation, RHMs assume that the measurement source is an element 
of a scaled version of the shape boundary, where the scaling factor is characterized by a particular 
probability distribution, i.e., the scaling factor is modeled as a one-dimensional random variable. 

Definition 2 (Random Hypersurface Model (RHM)). The measurement source z k t £ rn k + S(p ) 
on an extended object with star-convex shape <S(p,) C B 2 and center rn k £ <S(jO is generated 
according to an RHM if 

z k ,i ^m k + *k,i ■ s{p k ) , 

where S(p ) denotes the boundary of S(p ) and s kj i £ [0, 1] is a one-dimensional random variable. 

The restriction to star-convex shapes ensures that S(p ) = Usefo i]{ s ' ^(P k )}- ^ n ^ ac ^ scaling 
the object boundary corresponds to a straight-line homotopy from the object center to the object 
boundary. Note that all scaling factors s k ,i are mutually independent because measurements are 
assumed to be mutually independent. 

Definition [2] does not specify where the measurement source z k i lies on the scaled boundary. In 
general, it is possible to consider it as an unknown fixed parameter (functional model) or to assume 
it to be drawn from a probability distribution (structural model). These two models are also widely- 
used in curve fitting [32J. For the structural model, an RHM becomes a spatial distribution model 

mm- 



Remark 2. The implicit representation motivated in Section 3.1 is based on the functional model. 
Probabilistic information about the location of a measurement source on the boundary is not 
explicitely encoded in the implicit function. The shape boundary can be written as 

S(p k ) = {z \z G E 2 and g(z, p_ k ) = 0} , (4) 

and 

s k,i ■ S{p k ) = {z | z G R 2 and g(z,p_ k , s k ,i) = 0} (5) 

is the scaled shape boundary. 

3.3. Probability Distribution of the Scaling Factor 

The probability distribution of the scaling factor depends on the distribution of the measurement 
sources on the extended object. The following theorem says how the scaling factor is distributed 
in case the measurement sources are uniformly distributed on the surface of a two-dimensional 
extended object. 

Theorem 1. If the measurement source z£~R 2 is uniformly distributed over the star-convex region 
5cl 2 with center [0, 0] , the corresponding squared scaling factor is uniformly distributed on the 
interval [0, 1]. 

Proof. The cumulative distribution function F(s) of s turns out to be 

F(s) =p(s<s)=p(zEs-S)= A 2ea(,sf ) = g ' ( 6 ) 

for s G [0, 1]. Additionally, F(s) = for s < and F(s) = 1 for s > 1. The cumulative distribution 
of u := s 2 given by F(u) = p(s < y/u) is the cumulative distribution function of the uniform 
distribution on [0, 1]. □ 

According to Theorem [TJ it is reasonable choose a uniform distribution for the squared scaling 
factor. 

3.4- General Procedure for Extended Object Tracking with RHMs 

In order to derive a state estimator for an extended object based on an RHM, the following 
steps are to be performed: 

• A suitable shape and a shape parameterization have to be determined. 

• The implicit equation (I3| has to be formed. 

• A state estimator has to be derived based on the implicit measurement equation defined by 
(§ and (§. 

Based upon the above steps, we will derive particular Gaussian estimators for ellipses and free-form 
star-convex shapes in this article. 



4. Formal Gaussian State Estimator for Extended Objects 

In the following, the notation of a formal (Gaussian) Bayes filter for the state x k based on 
the previously discussed models is introduced. Particular implementations based on RHMs are 
presented in the next two sections. 

We denote the probability density for the parameter vector x k after the incorporation of all 
measurements up to time step k — 1 plus the measurements y ,.. . ,y with fi(x k ) for < I < 
n k . In this article, we focus on Gaussian state estimators so that all probability densities are 
approximated with Gaussians, i.e., fi{x k ) & M(x_ k — (j/P , Yi k t ), where n x is the mean and T, k t the 
covariance matrix. 

Time Update. The time update step predicts f n . k _ 1 (x k _ l ) to the next time step, i.e., it determines 
fo(x k ). As we focus on linear system models Q, the time update can be performed with the 
Kalman filter formulas, see for example pQ. 

Measurement Update. The prediction fo(x k ) is updated with the set of measurements {y, ,1^ 
according to Bayes' rule. Because the measurement generation process is assumed to be independent 
for consecutive measurements, they can be incorporated recursively according to 

fi(x k ) = a k ,i ■ f{y kl \x k ) ■ fi-i{x k ) for 1 < I < n k , 

where f(y,,\x k ) is a single measurement likelihood function and a ky i is a normalization factor. 

Note that the order of the measurements for a particular time step is irrelevant, because they are 
generated independent. However, the processing order may matter if approximations are performed 
(see Remark [3]) . 

5. Elliptic Shapes 

Elliptic shape approximations are highly relevant for real world applications as many targets are 
approximately elliptic. Even in case of high measurement noise and few available measurements, an 
elliptic shape approximation can be estimated. In this section, an implicit measurement equation 
is derived for elliptic shapes based on a RHM and subsequently a Gaussian state estimator is 
developed. 

5.1. Parameterization of an Ellipse 

An ellipse is determined by its center m k 6 R 2 and a positive semi-definite shape matrix 
Afc G R, 2x2 . Based on the Cholesky decomposition (A^.)" 1 = L^L^T, where 



a k 
c k b k 

a suitable vectorized parameterization p? 1 := [0^,6^,0^] of A^ can be defined. With x 
[jjlfc i (2±fc) T > (P?) ] j the implicit shape function becomes 

9 e \z kjl ,xt) := (z ktl - m k ) T ■ A k x ■ (z k>l - m k ) - 1 . 



(7) 



5.2. Implicite Measurement Equation 

Thanks to the chosen representation of an ellipse, the implicit function for the scaled version 
of <S(p?') with scaling factor s k j turns out to be 

9 cl Uk,hxf, a k ,i) ■= (z k)l - m k ) T • A^ 1 • (z Kl - m k ) - s 2 ktl . (8) 

Equation ([8]) and (TlJ) specify together an implicit measurement equation, which in this case coincides 
with the problem of fitting an ellipse to noisy data points plus the additional random scaling factor. 

5.3. Gaussian State Estimator 

Approaches based on the Kalman filter for state estimation with implicit measurement equations 
such as ([8]) and (TO) are well-known in literature. Typically, the implicit measurement equation is 
linearized around the measurement and state in order to render it explicit, see for example [33J. 
Here, we propose to perform algebraic reformulations followed by a statistical linearization around 
the measurement source as described in the following. When ph would be linear, we could rewrite 
the problem directly as an explicit measurement equation by plugging (fTl) into pi). As pi) is 
nonlinear, an exact reformulation is not possible. Nevertheless this reformulation can be performed 
approximately. For this purpose, the first step is to plug into ^, i.e., 

9 cl (y k j,x e k \s ki i) = g(z kil + v k}l ,p£,s kj i) . (9) 

After some minor simplifications of (19]) and exploiting pi), the following measurement equation 

h e \xt,v k>l ,s k , h y kl ) = (10) 

with a pseudo-measurement is obtained, where 

h c \x? k \v k>h s kh y kl ) := 9 cl (y_ kl ,P k ,s k ,i) - 2(z k>l - m fc ) T A^ V fc ,z - «*,j A fc 1 v k ,l ~ s l,i ( n ) 

maps the state x| , the measurement noise v k i, the scaling factor s kj i, and the measurement y 
to the pseudo- measurement. 

However, the unknown measurement source z k ; still occurs in (10). The basic idea here is 
to substitute z k t — m k in ( 10 ) with a proper point estimate. The easiest way to obtain a point 



estimate for z k l — rn k is to consider the ellipse specified by the mean of the previous estimate, i.e., 
\i? k l _ 1 and use the point with the smallest distance from the conic to the measurement y as a 
point estimate for obtaining a point estimat for z k t — m k . 



The measurement equation (10) can directly be used within Gaussian state estimators such 



as the Unsented Kalman Filter (UKF) |34j or analytic moment calculation [35], which essentially 



performs a statistical linearization of (10). More precisely, the joint density of the state and 
the predicted pseudo-measurement is approximated with a Gaussian. Then, the Kalman filtering 
p] formulas can be used for calculating the updated estimate, i.e., given the previous estimate 
fl-l(x k ) = M(x k - //^_ 1 ,S| )i _ 1 ), the updated estimate fi(x k ) = N(x k — //? ,,S^j) with the 
measurement y results from the Kalman filter equation: 

£ f| = ^_ 1 + Sf(Sf)- 1 (0-/xJ) , (12) 

Ex _ x^x x^xh/\r^hh\ — lx^hx /io\ 

k,l — ^kl-l ~ ^k \^k I ^k i \ 16 ) 



5 For sake of clarity, the index "el" is omitted for the mean and covariance matrices of the estimates. 



where is the predicted pseudo- measurement, T,f, is the covariance between the pseudo-measurement 
and the state, and S^ is the variance of the predicted pseudo- measurement. In this manner, a 
statistical linearization around the measurement source is performed. At this point, it is important 
to note that \xl and X£ do not depend on the unknown measurement source z k ; and, hence, the 
error made due to substituting it with a point estimate is rather negligible. 

Remark 3. Due to the nonlinear measurement model and the performed approximations, the order 
of the measurement processing matters. However, we observed that the differences arising from 
different processing orders is rather negligible. Besides, it is possible to perform a batch processing 



of all measurements by means of stacking the single measurement functions (10). By this means, 
the estimation accuracy can be slightly increased. 

Remark 4. Due to the statistical linearization of the implicit measurement equation, approximation 
errors are introduced. As a consequence, the resulting estimation quality highly depends on the 
parameterization of the ellipse and of the particular form of g el {y,,,x'k^ s k,i)- For example, we 
observed that multipliying both sides of (19]) with 1/ trace (L^L^ ) improves the estimation quality. 
Remark 5. Note that the actual likelihood function f(y, ,|a±|) can be derived from the measurement 



equation (17), however, it is not required explicitely when performing statistical linearization. 



6. Star-Convex Shapes 

When the measurement noise is low compared to the target extent, it may be possible to extract 
more detailed shape information than an ellipse from the measurements. For this purpose, an RHM 
for estimating and tracking the parameters of a star-convex shape approximation is presented in the 
following. A detailed target shape approximation is of high value for many higher-level problems 
such as target classification, track management, and sensor management. Last but not least, a 
more detailed shape estimate results in a better estimation quality for the kinematic state of the 
target as it is a more precise model of the reality. 

6.1. Parametrization of Star- Convex Shapes 

A star-convex shape S(p s ^) can be represented in parametric form with the help of a radius 
function r(p s £, </>), which gives the distance from the object center to a contour point depending on 
the angle <fi and a parameter vector p?° ( see Fig.[4J). A suitable (finite dimensional) parameterization 
p!° of a radius function is given by the first Nf Fourier coefficients that define a Fourier series 
expansion [361 EZ], i- e -i 



r (p-0) = !<4 O) + Yl (a^cosO^ + ^sin^)) = R(0) . p> 



where 



R(4>) = [^,cos((j)),sm((f)),. . . ,cos(Np(j}),sm(NF(j))] , and 

„sc _ r (0) (1) 6 (1) (N F ) b (N F ),T 

Ek ~ Ik ' a k ' °k ' ' ■ ■ a k i °k J 



Fourier coefficients with small indices encode information about the coarse features of the shape 
and Fourier coefficients with larger indices encode finer details. 
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(a) Star-convex shape. 

Figure 4: Representation of a star-convex shape with a radius function. 
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#.£?. Implicit Measurement Equation 

With ai| c := [m|\ (^) T , (p!°) T ] and the implicit representation of star-convex curves [36], we 
obtain 

5 SC Ufc,b^ C ) = ll?21-^|| 2 -^(p fa fc C ^(^-^)) 2 . ( 14 ) 

where Z(m — z k {) denotes the angle between m — z kl and the x-axis. The scaled version of the 
shape boundary is specified by 

9 s %z kp xf,s k: i) = \\m-z kl \\ 2 -sl r r(^ k c ,Z(rn-z kjl )) 2 , (15) 

Again, (15) specifies together with (fTl) an implicit measurement equation. 

6.3. Bayesian State Estimator 

A measurement equation can be derived by plugging into fl~5| ) 

g sc {y kl ,x k ,s kt i) = g{z Kl + v Kh ^,s Kl ) . (16) 

In order to avoid the treatment of uncertain angles, we propose to replace the angles in occurrences 
r(p? c , •) with a point estimate (p k> i, i.e., we assume r(p, c , Z(m — y )) w r(p sc , /.(m—(z k i+v kl )) ~ 

Remark 6. A proper point estimate is given by the most likely angle <j)k,i- In case of isotropic 
measurement noise, this point estimate (f> kt i is given by the angle between the vector from the current 
shape center estimate H%ri to the measurement y and the x-axis, i.e., 4>k,l '•= ^-(Vui ~ Mil\_i)- 



Based on the point estimate and rn k — z kl « s k j ■ R(cj>k,i) " P s k " §{ < t > k,i)-, where e( 



^WJ 



[cos^fc,/, sin^fe^] , (16) can be simplified to the following measurement equation 

h sc (x k ,Vk,h s k,i,y kl ) = 
with 



(17) 



h sc (xT,v k:l ,s k j,y kl ) :=<HI R 0M-P^II + 2s M R(</> M )p_;: c e(<M w*,j + IK 



ijfci 



m 4 



which maps the state a;? , the measurement noise ;y fc ;, the scaling factor s^/, and the measurement 
y to a pseudo- measurement 0. 
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(a) Target 1: Elliptic shape, (b) Target 2: Aircraft-like (c) Target 3: Group target. 

shape. 

Figure 5: Targets used for the evaluation. 



Remark 7. An alternative derivation of (17) based on the explicit representation of the shape with 
the radius function can be found in [351 . 



(x k ) = N{x k ~ Eki-i'^k.i-i)' tne posterior density fi(x k ) = M{x_ k - 

T, k i) having received the measurement y can be calculated with a Gaussian state estimator 

or analytic moment calculation [35J for a closed-form measurement update. 



For a given density //_i 
such as the UKF 13 



See also (12) and (13) in this context. Just as for ellipses, the order of the measurement processing 



matters (see the discussion in Remark k3|). 



7. Evaluation 

RHMs for elliptic and star-convex shapes are evaluated by means of both stationary and mov- 
ing extended objects r\ For both elliptic shapes and star-convex shapes the UKF |34] is used for 
performing the measurement update. For elliptic shapes, the squared scaling factor is modeled as 
a Gaussian distribution with mean 0.5 and variance 1/12 (i.e., the first two moments of a uniform 
distribution). For star-convex shapes, the scaling factor is modeled as Gaussian distribution with 
mean 0.7 and variance 0.06. 

1.1. Stationary Extended Target 

In the first scenario, we consider an extended object with fixed position and shape. From the 
target, 300 measurements are received sequentially, i.e., a single measurement per time step (n k = 
1). Simulations are performed with the three different target types depicted in Fig. pi where 
measurement sources are drawn uniformly from the target surface/group members. 

The parameters of the ellipse are a priori set to a Gaussian with mean [0.5,0.5, 1.6, 1.6,0.6] 
and covariance matrix diag(3, 3, 0.5,0.5,0.5), i.e., an uncertain circle with radius 1.2 and center 
[0.5,0.5] . For the star-convex shape approximation, 15 Fourier coefficients are used, where the 
shape parameters are a priori set to a Gaussian with mean [0.5, 0.5, 3, 0, . . . , 0] and covariance 
matrix diag(0.7, 0.7, 0.1, 0.1, . . . ,0.1), i.e., an uncertain circle with radius 1.5 and center [0.5, 0.5] . 



J Source code for RHMs is available at http://www.cloudrunner.eu. 
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(a) Low noise level. 
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(c) High noise level. 
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Figure 6: Simulation results: Example measurements for a particular run and point estimates for 
the shape averaged over 20 runs. Simulations are performed with low measurement noise level 
TV 1 = diag(0.6 2 ,0.6 2 ), medium measurement noise level Tit 1 = diag(l, 1), and high measurement 
noise level £j£ x = diag(1.4 2 , 1.4 2 ). 
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(a) Low noise level. 
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(b) Medium noise level. 
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(c) High noise level. 





a 




2 4 
x/m — > 



Figure 7: Simulation results: Example measurements for a particular run and point estimates 
for the shape averaged over 20 runs. Simulations are performed with low measurement noise 
level Ti" kl = diag(0.3 2 ,0.3 2 ), medium measurement noise level Ti" kl = diag(0.4 2 , 0.4 2 ), and high 
measurement noise level Y> v kl = diag(0.6 2 , 0.6 2 ). 
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(a) First measurement. 



(b) 10 measurements. 



(c) 30 measurements. 



(d) 50 measurements. 



Figure 8: Example run demonstrating the sequential incorporation of measurements (Noise level: 
S^ = diag(0.3,0.3)). 

The estimation results after the 300 measurements are depicted in Fig. [6] and Fig. [7]for different 
measurement noise levels. The shape estimates are averaged over 20 Monte-Carlo runs. In order to 
illustrate the magnitude of the measurement noise, the measurements of a particular run are also 
plotted. It is important to note that this is just done for visualization as the estimator incorporates 
the measurements recursively. Fig. [8] depicts a single example run for star-convex shapes in order 
to show the evolution of the shape estimates with an increasing number of measurements. 

1.2. Moving Extended Object 

In the second scenario, the aircraft-shaped target shown in Fig. [9a| moves along the trajectory 



depicted in Fig. 9b The measurement sources are drawn uniformly from the target surface. The 



magnitude of the measurement noise varies from measurement to measurement in order to simulate 
different sensors or different target-to-sensor geometries. The covariance matrix of the measurement 
noise is Ej^ = diag(0.2 2 ,0.2 2 ) with probability 0.75 and Ej^ = diag(0.4 2 , 0.4 2 ) with probability 
0.25. The number of measurements received per time instant is given by n^ = n* k + 1, where n* k is 
a Poisson distributed random variable with mean 4 for the ellipse and mean 7 for the star-convex 
shape. 

We employ a constant velocity model for the temporal evolution of the target center [T] and 
a random walk model for the shape parameters. Hence, the state vector to be tracked is x k = 
[m k , m k , (p s , c ) T ] , where m k is the center, m k is the velocity vector, and p, c are the shape 
parameters. As the center is assumed to evolve according to a constant velocity model, the system 

I 2 TI 2 1 



matrix in ([2]) is A^ 
with dimension d - 



diag(A™, 1^), where A 



with T = 1 and 1^ is the identity matrix 



I 2 

5 for ellipses and d = 11 for star-convex shapes. The system noise is zero- mean 

T3 I 2 

0.0015 and q<i = 0.005. For the star-convex shape q± 



Gaussian noise with covariance matrix C k 



ellipse, we have q\ 



diag(gi -I* C£>) withCJ" = «2 



0.0001 and q 2 



For the 



0.003. 



The estimated shapes (averaged over 20 time steps) are depicted in Fig. 10 and Fig. 11 for 



two snippets of the trajectory. The results show that the shape of the extended object is tracked 
precisely, even when the shape changes its orientation. 
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Figure 9: Extended object and its trajectory. 
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(c) Example measurements. 



(d) Shape estimates. 



Figure 10: Tracking an extended object with an RHM for ellipses: Example measurements from a 
particular run and average shape estimates. 



Note that in the simulations with elliptic shapes the number of measurements per time step is 
lower than for star-convex shapes. A star-convex shape approximation can only be extracted well 
in case the measurements carry enough information, i.e., enough measurements with rather low 
noise are available. 



The example measurements in Fig. 10 and Fig. 11 also emphasize that naive approaches for 
estimating a shape would be bound to fail, e.g., directly computing an enclosing shape of the 
measurements is infeasible because the measurements are noisy and only a couple of measurements 
are available per time step. 



8. Conclusions and Future Work 

This article considered the problem of estimating a shape approximation of an extended object, 
which gives rise to several measurements from different spatially distributed measurement sources. 
For this purpose, a novel approach for modeling extended objects called Random Hypersurface 
Model (RHM) was introduced that allows to derive a functional relationship between the measure- 
ments and shape parameters. We presented particular RHMs for elliptic and free-form star-convex 
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(a) Example measurements. (b) Shape estimates. 



(c) Example measurements. 



(d) Shape estimates. 



Figure 11: Tracking an extended object with an RHM for star-convex shapes: Example measure- 
ments from a particular run and average shape estimates. 

shapes and derived measurement equations for which standard Gaussian state estimators can be 
used. Nevertheless, estimating a detailed star-convex shape approximation is only possible when 
the measurement noise is rather low compared to target extent and enough measurements per time 
step are available. If this not the case, a basic shape such as an ellipse is more suitable. Hence, a 
mechanism for adapting the complexity of the used shape description is desired. The capability of 
estimating free- form shape approximations paves the way for new applications, e.g., classification 
based on the shape and group splitting detection. 
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